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It is known that in a water solution with multivalent counterions (Z-ions) two likely charged 
macroions can attract each other due to correlations of Z-ions adsorbed on their surfaces. This 
"correlation" attraction is short-ranged and decays exponentially with increasing distance between 
macroions at characteristic distance A/2n, where A is the average distance between Z-ions on the 
surfaces of macroions. In this work, we show that an additional long range "polarization" attraction 
exists when the bare surface charge densities of the two macroions have the same sign, but are 
different in absolute values. The key idea is that with adsorbed Z-ions, two insulating macroions can 
be considered as conductors with fixed but different electric potentials. Each potential is determined 
by the difference between the entropic bulk chemical potential of a Z-ion and its correlation chemical 
potential at the surface of the macroion determined by its bare surface charge density. When the 
two macroions are close enough, they get polarized in such a way that their adjacent spots form a 
charged capacitor, which leads to attraction. In a salt free solution this polarization attractive force 
is long ranged: it decays as a power of the distance between the surfaces of two macroions, d. The 
polarization force decays slower than the van der Waals attraction and therefore is much larger than 
it in a large range of distances. In the presence of large amount of monovalent salt, the polarization 
attraction decays exponentially at d larger than the Debye-Hiickel screening radius r^. Still, when 
A/2n <C d <C Ts, this force is much stronger than the van der Waals attraction and the correlation 
attraction mentioned above. The recent atomic force experiment has shown some evidence for this 
polarization attraction. 

PACS numbers: 61.20.Qg,61.25.Hq,82.45.Gj 



I. INTRODUCTION 

Water solutions of strongly charged colloidal parti- 
cles (macroions) vi^ith multivalent [Z valent) counteri- 
ons (Z-ions) are important in physics, chemistry, biology 
and chemical engineering. Colloidal particles, charged 
Hpid membranes, DNA, actin, and even cells and viruses 
are examples of different macroions. Multivalent metal- 
lic ions, dendrimers, charged micelles, short DNA and 
other relatively short polyelectrolytes like spermine can 
play the role of Z-ions. Several interesting, counterintu- 
itive phenomena have been discovered in such a system, 
such as charge inversion, which have attracted signifi- 
cant attention (see review paper Ref. Q and references 
therein). Charge inversion happens in a water solution 
when a macroion binds so many Z-ions that its net charge 
changes sign. This phenomenon is especially important 
for gene therapy. To deliver DNA into the cell, charge 
of negative DNA (playing the role of macroion) should 
be inverted by positive Z-ions to approach a negatively 
charged cell membrane. Here we actually assume that the 
membrane is weakly charged and therefore its charge is 
not inverted. The question is what happens if the charge 
of the membrane is also inverted by positive Z-ions. Can 
DNA still be attracted to the membrane? 

A similar question was put forward by the recent 
atomic force experiment ,2,] designed to verify the the- 
ory of charge inversion based on strong correlations be- 
tween Z-ions Qjll'l In the experiment, forces between 
a negatively charged spherical macroion attached to the 
cantilever (probe) and a positively charged surface were 
measured at different concentrations of positive Z-ions 



(see Fig. ^) 0j . At small concentrations of Z-ions, the 
probe is attracted to the surface. With increasing con- 
centration of Z-ions, the charge of the probe gets inverted 
by Z-ions, and the measured force at large d becomes re- 
pulsive, where d is the distance of closest proximity from 
the probe to the surface. The critical concentration of 
Z-ions, where this happens is in reasonable agreement 
with the prediction of Ref. However, an interesting 
new feature of the repulsive force was observed. With 
decreasing d, the repulsive force reaches a maximum at 
d ~ 100 A (which is roughly equal to the Debye-Hiickel 
screening radius of the solution) and start to decrease. 
This suggests the existence of a competing attraction. 

One can also consider a different experimental setup 
in which the probe and the surface are likely charged and 
both adsorb Z-ions (Fig.^). When the concentration of 
Z-ions is high enough so that the charges of the surface 
and the probe are both inverted at large d, it is inter- 
esting to find out whether the force may be attractive. 
Preliminary experiments showed such attraction ^ . No- 
tice that this setup brings us back to the original question 
of attraction of charge-inverted DNA to charge-inverted 
cell membrane. 

Motivated by these questions, in this paper we study 
the interaction between two different macroions in the 
presence of a large concentration of Z-ions. In the main 
part of the paper, we focus on the case of two bare nega- 
tively charged macroions and positive Z-ions correspond- 
ing to Fig. ^p, since it is pedagogically simple. Only in 
Sec. we discuss the case of oppositely bare charged 
macroions (Fig.QJt) in the connection with Ref. 0- We 
assume that the valence of Z-ions, Z ^ I, but still many 
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FIG. 1; Schematic illustration of two experimental setups for 
atomic force measurements of likely charged macroions. The 
force between the probe (the big sphere with negative charge) 
and the surface is measured at different concentrations of Z- 
ions (small spheres with positive charges), (a) In the setup 
of Ref. 0, the bare surface and the bare probe are the op- 
positely charged so that positive Z-ions are only adsorbed to 
the probe and invert its charge, (b) In the setup of Ref. Q, 
the bare surface and the bare probe are likely charged so that 
positive Z-ions are adsorbed to both of them, again making 
them likely charged at large d. 

Z-ions are needed to neutralize one macroion. The two 
macroions are spherical but different in their bare sur- 
face charge densities. We will see that this is crucial to 
produce the attraction. Actually this is the reason why 
this attraction was not reported before, when the focus 
was on identical macroions P]. In calculation, we focus 
on the two limiting cases, R\ = i?2 and R\ <^Ri, where 
R\ and Ri are the radii of the two macroions. 




FIG. 2: Z-ions are strongly correlated on the surface of the 
macroion. Their short range order is similar to that of a 
Wigner crystal with lattice constant A. 

Before we discuss the mechanism of this attraction, 
let us first briefly review the theory of charge inver- 
sion 0. Let us consider a water solution with one neg- 
atively charged macroion and many positively charged 
Z-ions. Due to the Coulomb interaction, Z-ions are ad- 
sorbed to the surface of the macroion. On the surface, 
they strongly repel each other with energy much larger 
than koT and form a two dimensional strongly corre- 
lated liquid with short range order similar to a Wigner 
crystal (WC) (Fig. [21). When a new Z-ion approaches the 
macroion, it repels already adsorbed Z-ions and creates 
a negative spot for itself. One can view it as an elec- 
trostatic image of the Z-ion, similar to the image on a 
conventional conducting surface. Attraction to the image 
leads to an additional negative chemical potential /ic (c 
stands for "correlation" ) for each Z-ion on the surface of 



the macroion. As a result, when the concentration of Z- 
ions in the solution is large enough, so many Z-ions are 
adsorbed that the net charge of the macroion is inverted 
and becomes positive Q. 

As pointed out in Ref. for an already neutral 
macroion, adsorption of additional Z-ions can be viewed 
as charging a conductor of certain capacitance with a 
fixed electric potential 0. This potential is determined by 
the difference between the chemical potential of Z-ions in 
the solution and the correlation chemical potential /ic on 
the surface. It is a constant all over the macroion surface 
because the net surface charge density of the macroion 
changes very little and therefore /ic is a constant (see 
Sec. m for detail discussion) . Actually, this electric po- 
tential is similar to the surface potential of a colloidal 
particle determined by its "potential-determining ions" 
given by the Nernst equation The only difference is 
that in the present case, [i^ originates from correlations 
of Z-ions instead of chemical bonding. 

The model of conductors is convenient to discuss the 
interaction and understand the possible attraction be- 
tween two different macroions in the presence of Z-ions. 
Actually, as we know from electrostatics, if two conduc- 
tors are charged with different potentials, they may at- 
tract each other, even though the two potentials have 
the same sign. For example, consider two conducting 
spheres 1 and 2 with same size and potentials and (^i 
respectively and (/ii > 02 > 0. When the two spheres 
are far away, they repel each other. When we bring 
them closer, adjacent spots of two spheres form a charged 
capacitor with potential difference (\)\ ~ <i>2- In other 
words, spheres get polarized. The polarization attraction 
competes with the overall repulsion between two spheres 
(since qi > q2 > 0) and dominates at small distances be- 
tween spheres. Similar polarization leads to attraction of 
a conducting sphere to a close conducting plate holding 
at different potentials. 

Using the model of conductors, we find that as long as 
the bare surface charge densities of the two macroions are 
different, they always attract each other at small enough 
distances, even when their bare charges or net charges 
are of the same sign. This polarization attraction is even 
more appreciable in the case that the sizes of the two 
macroions are very different as in Fig. ^ We also find 
that the attractive polarization force Fp decays with in- 
creasing distance as a power law. In the case of Fig. ^ 
at d <C r-s , we get 

F,^-^(0r-0.f, (1) 

where is the Debye-Hiickel screening radius, D = 80 
is the dielectric constant of water and Ri is the radius of 
the spherical probe. 

We emphasize that the polarization attraction dis- 
cussed here is different from two standard attractive 
forces. The first one is the well known van der Waals 
attraction used in standard DLVO theory . In the case 
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of Fig. n it is given by 
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where H is the Hamaker constant. It is clear that Fy^w 
decays with d faster than Fp. The ratio of these two 
forces is 
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For the typical if = 1.0 x 10'^° J = 2AkBT and reason- 
able — (^2! — OAkBT/e (see Sec.mi, this ratio is larger 
than unity if c? > 6 A. 

The second competing force is the short range attrac- 
tion between two macroions due to correlations of Z-ions 
on their surfaces |^. In the spot where two macroions 
touch each other, the surface density of Z-ions is doubled 
and the correlation energy is gained. This correlation 
attraction is related to Wigner-crystal-like arrangement 
of Z-ions on the surfaces and therefore decays exponen- 
tially with increasing distance as e'^'^'^M^ where A is the 
"lattice constant" of the Wigner crystal (see Fig.|21). For 
d ^ A/2TT, it becomes much weaker than the polarization 
attraction given by Eq. ^ (we assume that ^ A/2tt). 
In this sense the polarization attraction is a long range 
force. 

One can understand the polarization attraction from 
another point of view, i.e., from the concept of contact 
electrification. As well known, when two different solids 
contact to each other in vacuum, due to the difference in 
their work functions, certain amount of electrons move 
from one material to the other, developing a electric volt- 
age which stops further charging. This contact electrifi- 
cation leads to the well known Coulomb attraction be- 
tween them 0. Here a contact between two objects is 
necessary to eliminate the kinetic barrier for electrons 
and makes electrification possible during time of exper- 
iments. If we wait long enough, the electrification can 
happen even without a close contact. In our case of 
macroions in water, Z-ions play the role of electrons and 
the absolute value of chemical potential fic of a Z-ion on 
the surface of the macroion plays the role of the work 
function. To keep the same electro-chemical potential of 
Z-ions, since /Xc is different for two macroions, the electric 
potentials must also be different, which leads to attrac- 
tion. The difference in our system is that the kinetic 
barrier for Z-ions is relatively small and the equilibrium 
of Z-ions is easily achieved through the solution. 

Similar electrification phenomenon has been studied 
m Ref. ^ in a toy model of two negative spheri- 
cal macroions exactly neutralized by Z-ions. The two 
spheres had the same radius, but substantially different 
bare charges. It was shown that under conditions of total 
neutrality one sphere becomes undercharged (negative) 
and the other is overcharged (positive) and, therefore, 
they attract each other. 

In the present paper, we discuss a more generic and 
realistic situation: there is certain concentration of Z- 
ions in the solution (a fixed chemical potential of Z-ions) 



and the macroions can be both either undercharged or 
overcharged by them. In spite of the same sign of their 
net charges these macroions attract each other at small 
distances because of local polarization around the points 
of closest proximity. 

In the presence of large amount of monovalent salt, 
the Coulomb interaction is effectively truncated at the 
Debye-Hiickel screening radius rs- When = (j)2, the 
two macroions with adsorbed Z-ions repel each other as 
described by the standard DLVO theory f5l|. However, 
when 01 ^ (j)2, i.e., when the two macroions are differ- 
ent, the polarization attraction appears. When t/fi and 
(f>2 have the same sign, it decays as e"^'^/''^ at distances 
larger than Vg- Still, for A/tt <C d ^ r^, this attraction is 
much stronger than both the van der Waals attraction Q 
and the correlation attraction P|, which is proportional 
to e-2'^''/'4. 

This paper is organized as follows. In Sec. ^ we de- 
scribe the interaction between two different macroions in 
the presence of Z-ions and show that it is similar to the 
interaction between two conductors with different poten- 
tials. In Sec. mil we discuss the interaction of two con- 
ducting spheres and derive the power law of the attractive 
force. In Sec. IIVI we take into account of the effect of 
screening by monovalent salt. In Sec. we generalize 
our theory to the case when only one macroion adsorbs 
Z-ions in the connection with the experiment in Ref. 
(see Fig.^). We conclude in Sec. IVTl 



II. THE MODEL OF CONDUCTORS FOR TWO 
MACROIONS IN THE PRESENCE OF Z-IONS 



In this section we show that the interaction between 
two macroions in the presence of Z-ions can be consid- 
ered as the interaction between conductors with different 
electric potentials. Let us first consider a single spheri- 
cal macroion with radius R and bare charge — Q < in a 
water solution with Z-ion concentration n. The bare sur- 
face charge density of the macroion is —cr = —Q/AttR^. 
As we discussed in Sec. 13 Z-ions on the surface of the 
macroion strongly repel each other to form a strong cor- 
related liquid, similar to a structure of Wigner crystal 
in the short range (see Fig. The chemical potential 
related to this correlation is dominated by its low tem- 
perature expression, which can be estimated from the 
Coulomb interaction inside each WC cell (see Eq. (9) in 
Ref. '3| for the full finite temperature expression of ^.c) 
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Dr D 
Here s is the surface density of Z-ions on the macroion. 



and r = A\J \/3/27r is the radius of the WC cell (see 
Fig. EJ). It satisfies irr'^s = 1. 

The equilibrium condition of a Z-ion in the solution 
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then reads 



ZeQ* 
DR 
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(5) 



where Q* is the net charge of the macroion combined 
with adsorbed Z-ions. The left hand side is the electro- 
static energy of a Z-ion at the surface of the macroion 
plus the correlation chemical potential. The right hand 
side is the difference between the entropic parts of the 
chemical potentials of a Z-ion in the solution and at the 
surface. As we will see below, this equilibrium condition 
is the key for the analogy with conductors. 

When n increases, the number of Z-ions adsorbed to 
the macroion also increases. At some critical n — uq, Q* 
becomes zero. According to Eq. we have 
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Here sq — a/Ze is the surface density of Z-ions to 
neutralize the macroion. The corresponding radius of 
WC cell is tq. In the case we are interested, /i(so) ~ 
Z^e^/Dro 3> fc^T (this defines "strongly charged" Z- 
ions). Therefore uq is an exponentially small concentra- 
tion which is easy to reach in experiments. For n > uq, 
more Z-ions come to the macroion and Q* becomes posi- 
tive. According to Eq. ((SJ, when n is so large that the en- 
tropy term is completely negligible, the maximum value 
of inverted Q*, Q^axi achieved Q, 



Q:„,, = 0.8 v/gZe. 
It is calculated by assuming 



(7) 



(8) 



This assumption is self-consistent, since for Q ^ Ze we 

have {Sraax " So)/sO = Q*max/Q = 0-8y/Ze/Q < 1. 

Now we can clearly see the similarity between a neu- 
tralized macroion and a neutral conductor. For a given 
n > no, the charging process from Q* = to Q* > can 
be viewed as charging a conductor with fixed potential. 
Indeed, the equilibrium condition Eq. (jSJ can be written 
as 



with 



Q* = DR(t>, 



\tJ-c{so)\ - ksTln^sl^^ /n) 
Ze 



(9) 



(10) 



Notice that (j) can be expressed through sq because of 
Eq. (|SJ). Eq. Q is exactly the same as the expression 
of charge for a conductor with capacitance C — DR and 
charged at fixed potential ip. For n < no, by similar dis- 
cussion, we have a conductor charged at < 0, providing 
that |Q*(n)| <C Q so that Eq. (|HJl is stiU satisfied. 

Now let us consider two macroions in a water solution 
with Z-ions. Clearly, when d is very large, if n is close to 



no for both macroions, the analogy to conductors holds. 
The two macroions have electric potentials ipi and (j)2. 
When d decreases, due to the interaction between two 
macroions, s changes and becomes nonuniform on their 
surfaces. When d is very small, s may become substan- 
tially different from sq at closest spots of the two surfaces. 
Then Eq. (|SJ| is not valid and 0i^2 start to change with 
d .IQj . In this paper, we stay in the limit of large sq so 
that |s — Sol ^ So- Accordingly, Eq. (jS)) is always valid 
and our approximation of fixed potentials is fine. 
Rewriting Eq. (|10|l for two macroions, we have 



Z^ 



fl,2 



(11) 



Here the subindexes 1, 2 represent two macroions respec- 
tively. si_2 are the number densities of Z-ions which neu- 
tralize two macroions, i.e., si_2 = 0'1,2/^e. Notice that 
01,2 are completely determined by n and si,2, i.e., the 
concentration of Z-ions and the bare surface charge den- 
sities of the two macroions. For a given solution of Z-ion 
concentration n, if ci 7^ (72, 4>i 1^ 4'2- For example, for 
typical values Z = +3, n = 1 mM, — cti — —0.75 e/nm^ 
and —(72 = —0.45 e/nm^, we have 4>i = OJlksT/e, 
02 = O.ST/csT/e and |0i - 021 = 0.40fcBT/e. Here 
Mc(si) = — 7.6fcsT and iJic{s2) — —h.lkBT are calculated 
using the full expression of fic given by Eq. (9) in Ref. 0]. 

In order to calculate the force between two macroions, 
we can imagine that they are conductors in ion-free water 
with potentials 0i and 02 supported by two batteries. 
Such conductors are not in equilibrium with each other 
while our macroions are in equilibrium. This however 
does not matter for the calculation of the force which is 
the same in both cases. Indeed, the force depends only 
on potentials and capacitance matrix of the system. 



III. ATTRACTION OF TWO CONDUCTING 
SPHERES WITH DIFFERENT POTENTIALS 

In this section, let us focus on the interaction of two 
conducting spheres 1 and 2 with radii i?i and R2 and 
potentials 0i and 02 in ion-free water (see Fig. As 
we explained in the previous section, this interaction is 
equivalent to the interaction between macroions covered 
by Z-ions. We define the the distance of closest proximity 
between the two spheres as d and the distance between 
the centers of the two spheres as L = d + Ri + R2 (see 
Fig. O. Below we use L or d alternatively for conve- 
nience. We focus on two limiting cases when i?i = R2 
and Ri <C i?2 , but our approach is applicable for any two 
spheres. 

We start from the total free energy of two conductors 
with fixed potentials 7j 

U{L) = -^Cn^l - ic220^ - Ci20i02. (12) 

Here Cn, C22, C12 are the self and mutual capacitances of 
the two spheres depending on L. Notice that in the case 
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FIG. 3: The interaction between two spherical macroions in 
the presence of Z-ions can be considered as the interaction 
between two conducting spheres with fixed but different elec- 
tric potentials. We focus on two cases: (a) The size of the 
two spheres are equal, (b) The size of sphere 2 is much larger 
than sphere f (the radius of sphere 2, R2, is not shown). For 
d <^ R2, the interaction between two spheres can be con- 
sidered as the interaction between sphere f with its "image 
sphere" inside sphere 2. 



of fixed potentials, the work done by the environment to 
charge the two conductors should be included in the total 
energy and therefore each term in Eq. (|12|l has a negative 
sign. The generic formula of the force between the two 
conductors is given by 



dU dCii 



dC 



22 



dC 



12 



dL dL 2 



dL 2 



dL 



M2. (13) 



As will be shown, the force is always attractive when 
4>i 7^ 02 and the two spheres are close enough. And it 
decays as a power law with increasing distance between 
spheres. 



A. Two spheres of the same size 

Let us start from a simple case when the radii of the 
two spheres are equal, — R2 ~ R (Fig. In this 
case, Cii = C22 and Eq. ifT^ is simphfied to 



F{L) 



dCu 



dL 



dCi2 
dL 



(14) 



According to the standard method of image charges, Cn 
and C12 can be calculated by considering a infinite series 



of images induced in the two spheres [9| . We have 



Cii = i:>i?sinhtV — — -, 

sinh(2n- l)i' 

Cr2 = -^i^sinh.^:-^. 
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Here a„(t) is defined as 

o-n (t) — n sinh t cosh nt — cosh t sinh nt . 



(16a) 
(16b) 

(17) 



Since a„(i) > for t > 0, we have dCn/dL < 

0, dCu/dL > 0. Also one can easily prove that 
\dCii/dL\ < \dCi2/dL\. These inequalities lead to some 
simple results. Firstly, when the signs of (jji and 02 are 
opposite, F < 0, i.e., the force is attractive, as it should 
be. Secondly, when (f>i — 01 (j)2 — 0, we have F < 0. 
Physically, the sphere with nonzero potential induces op- 
posite charge in the sphere with zero potential, therefore 
they attract each other. Finally, when (j>i = (j>2, F > 0, 

1. e., the force is repulsive. Indeed, the two spheres with 
same potential can be considered as a capacitor as a 
whole. Obviously, the closer the two spheres are, the 
smaller the capacitance is, and the higher the energy is 
(remember the negative sign in Eq. H12|l ). 

We are specifically interested in the case where the 
signs of 01 and 02 are the same but the magnitudes are 
different. Without losing generality, we consider the case 
of 01 > 02 > 0. Introducing a = 02/0i, we have 



FiL) 



1) 



dCii/dL 
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dCi2 



The critical value of a at which F 



dL 

is therefore 



02. (18) 
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(19) 



When a < ckc, -F < 0, and vice versa. In Fig. 01 we plot 
Q!c as a function of L/2R by taking the first 1000 terms 
in the series of Eq. (|16|l . We see that the attraction is 
possible at any distance as far as the ratio of the two 
potentials is smaller than ac- The closer the two spheres 
are, the larger ac is, and the easier the attraction can be 
developed. 

In order to gain more intuition about the interaction, 
also determine the dependence of the attractive force on 
the distance analytically, we study asymptotic behaviors 
of the force in the case of 0i > 02 > 0. In the limit of 
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FIG. 4: The critical ac — cf>2/4>\, at which the interactive 
force between two spheres is zero, is plotted as a function of 
L/2R numerically. The minus (plus) sign under (above) the 
curve means that the interaction between the two spheres is 
attractive (repulsive) . 



L ^ 2i?, keeping the leading order terms in Eq. \ib\ . we 

get 



Cii = DR 



Ci2 — 



DR"" 



and 
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(21) 



The physical meaning is clear. When the two spheres are 
far away from each other, they interact like two point 
charges. In the zero order, the magnitude of the charge is 
determined by their own capacitance and potential, i.e., 
DR<pi and DR(f)2- This gives the first term in Eq. |(5lJ 
(one D in charges is cancelled out since the interaction 
itself is proportional to l/D). In the next order, each 
point charge induces opposite charge in the other sphere 
with the magnitude —DR^(f)/L following the standard 
method of images ^ . The interaction between the charge 
and its image gives the second term in Eq. H21|l . From 
Eq. 1)2 one can easily get ac — R/L. Therefore the 
attraction dominates only if (/)2 is much smaller than (pi. 
The attractive force is proportional to (R/L)^. 

The more interesting limit is when the two spheres are 
very close to each other, i.e., d R. In this limit, we 
have 



Cl2 



DR , 16(3^R 
m , 

4 d ' 

DR^ I3^R 



(22a) 
(22b) 



where /3 — 1.78 (see Appendix for the derivation of 
Eq. According to Eq. (df, we have 



Fid) 
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(23) 



Remarkably, in this limit Oc = 1 and the force is always 
attractive (the only exception is 4>i = 02 when the lead- 
ing order term given by Eq. (|^ is zero and one has to go 
to the next order). Therefore we conclude that for any 
small — 02 , two conducting spheres always attract each 
other at small enough distance. This originates from the 
fact that in the limit of d — > 0, Cn and C12 both diverge 
and diverge in the same way (see Ea. (|22|l ). 

It should be mentioned that for the interaction be- 
tween two macroions, the model of conductors is valid 
only for d ^ A. Indeed, instead of smeared charge dis- 
tribution on a real conductor surface, on the surface of 
the macroion, Z-ions are discrete and form a WC-like 
liquid with "lattice constant" A (see Fig. |21l. When d 
becomes comparable with A, the electric potential <j> gets 
a periodic component along the surface, which does not 
exist for a real conductor. In this paper, we always fo- 
cus on the case when d ^ A. Still, for Ze Q and, 
therefore, A R, there is a big window of A ^ c? <C i? 
in which the model of conductors works and Eq. H23|) is 
applicable. 



B. A small sphere close to a big sphere 

Now let us consider the interaction between a small 
sphere with radius Ri and a big sphere with radius i?2. 
We are interested in the hmit i?i ^ R2 (Fig. Eh)- We 
can solve this problem by a similar procedure as in the 
last subsection, starting from a formula like Eq. H15() but 
with both Ri and R2 in it 0. Instead of doing this 
complicated calculation, let us use physics intuition and 
look at the interesting limit d <C i?2- In this limit, the 
size of the large sphere, R2, becomes irrelevant to the 
problem, and the interaction can be considered as the 
interaction between a sphere with a metallic semi-space 
(see Fig.Et)). 

Let us first consider the case when 02 = 0. As well 
known in electrostatics, a point charge induces an image 
charge in a metallice semi-space which describes the in- 
teraction between the charge and the metal. Similarly, 
using a series of image charges, one can show that the in- 
teraction between a sphere and a metallic semi-sphere is 
equivalent to the interaction between the sphere and its 
"image sphere" induced in the metal (see Fig. ISb). The 
image sphere and the original sphere have the same size, 
and their positions are symmetric about the boundary 
of the metal. If the original sphere has fixed potential 

01, its image sphere has fixed potential — 0i. Since these 
two spheres are oppositely charged, we immediately see 
that a conducting sphere is always attracted to a metallic 
semi-space. Therefore spheres 1 and 2 attract each other. 

Now we consider the case when 02 ^ 0. This potential 
is equivalent to a charge sitting at the center of sphere 

2, with magnitude -Di?202- It induces an image charge 
in sphere 1. In the limit of i?2 ^ 00, this image charge 
is located at the center of sphere 1 with the magnitude 
— Z3i?i02 (see discussion of method of images in Ref. P|). 
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It is equivalent to add an potential — (/)2 to sphere 1. As a 
result, the potential of sphere 1 is renormalized to 4>i — 4'2- 
Therefore, in this particular case, the interaction between 
two spheres with potentials </>! and (j)2 is the same as 
two spheres with potentials — (j)2 and 0. Clearly, the 
attractive nature of the interaction between sphere 1 and 
its image sphere is not changed by the nonzero 02 (except 
the special case (pi = 02 when there is no interaction 
between the two spheres). 

Having established the attractive nature of the force, 
we now calculate it quantitatively. Since the interaction 
is essentially between sphere 1 and its image sphere, Cn, 
Ci2 and C22 can be expressed through a linear combi- 
nation of Caa and Cah, where Caa and Cab are self and 
mutual capacitances of sphere 1 and its image sphere 
given by Eq. \ib\ . We have 

Cii{d) - C„a(2d) - Ca6(2d), (24a) 

CM = Cab{2d)~Caa{2d), (24b) 
C22{d) = DR2+Caa{2d)~Cab{2d), (24c) 

Here we put the d dependence in each C to remind us that 
the distance between spheres 1 and 2 is d but the distance 
between sphere 1 and its image sphere is 2d (see Fig-Oa). 
Noticing that Cn — — C12 and dCn/dd — dC22/dd, 
using Eq. ((T!^ . we have 



Fid) 



dCnjd) (01-02)' 
dd 2 

_ rdCabi2d) dCaa{2d) 



dd 



dd 



.(25) 



Since dCaa/dd < and dCab/dd > as we discussed in 
the last subsection, the force is attractive for any given 
01 and 02 (except the special case 0i = 02 when there is 
no interaction between the two spheres). 

The expressions of forces are particularly simple in cer- 
tain limits. When Ri <C (i <C i?2, using Eq. (^0)) . we have 



Fid) 



DRl 
4d2 ' 



When d <C using Eq. (|^ . we have 



Fid) 



DRl 



if 



(26) 



(27) 



This result can be compared with Eq. H23|l if i?i — R. We 
see that the force is stronger by a factor of 2 in the case 
of two spheres with different sizes. This is simply due to 
the fact that the surfaces of the two spheres are "closer" 
and the capacitances diverge faster in the present case. 



IV. THE EFFECT OF SCREENING BY 
MONOVALENT SALT 

In a water solution of macroions and Z-ions, nor- 
mally there is also certain amount of monovalent salt. 



In most cases, the effect of screening by monovalent salt 
can be described by the linear Debye-Hiickel screening 
radius, (we discuss the possibility and effect of non- 
linear screening at the end of this section). What we 
discussed in the last two sections corresponds to the case 
rs ^ i?i , i?2 7 d. Now we would like to consider the oppo- 
site limit, A ^ <C i?2. We show that even though 
attraction is suppressed and loses to repulsion at d ^ , 
the two spheres still attract each other at d < . 



The method of images in the presence of 
monovalent salt 



Let us first discuss how the method of images is mod- 
ified in the presence of monovalent salt. We start from 
the simplest situation in which there is a point charge q 
close to a metallic semi-space in a water solution with 
monovalent salt described by screening radius r^. When 

— > CX3, it is well known that the equal potential con- 
dition at the boundary of the metal can be satisfied by 
putting an image charge q' = —q inside the metal at the 
position symmetric to q (Here and below, as a premise, 
water should also be introduced with the image charge 
to produce the same dielectric constant). In the pres- 
ence of monovalent salt, one can easily check that the 
same boundary condition can be satisfied by introducing 
the same image charge q' at the same position, providing 
that a virtual cloud of monovalent ions with the screening 
radius is produced together with the image charge. 





(b) 



FIG. 5: Images in a conducting sphere in the presence of 
monovalent salt with screening radius r^: (a) a point charge 
q induces an image charge q' — —Rq/L with a virtual cloud 
of monovalent ions with screening radius r's — Rvs/L; (b) a 
conducting sphere with potential is equivalent to a point 
charge q" = DRcjje^^^" with a virtual cloud of monovalent 
ions with screening radius r^. 

As a more complicated situation, let us consider a point 
charge q and a grounded conducting sphere in a water 
solution with monovalent salt described by (Fig. |3^). 
As well known ?|, when Ts ^ 00, the magnitude and 
position of the image charge are given by (see Fig. |31 for 
definition of all lengthes) 



R 



L' 



R^ 
L ' 



(28) 
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When Ts is finite, similar to the case of metalhc semi- 
space, a virtual cloud of monovalent ions should be cre- 
ated together with the image charge. Due to the special 
geometry of a sphere, as one can check, the screening 
radius of this virtual cloud is given by 



different from in the solution. 

Another relevant situation is a conducting sphere with 
fixed potential in a water solution (see Fig. ^p) . In this 
case, the potential in the solution can be solved exactly 
using linearized Poisson-Boltzmann equation and inter- 
preted by the method of images. Actually it is equivalent 
to the potential produced by a image charge q" at the 
center of the sphere. We can create a virtual cloud of 
monovalent salt together with the image charge as our 
convenience. The magnitude of q" is then determined 
not only by the boundary condition but also the virtual 
cloud we choose. As we will see in the next subsection, 
it is convenient to introduce inside the sphere and get 



from the leading order image charges only. To see this, 
one just remember that the image charges really repre- 
sent surface charge densities induced on the two spheres. 
Each image in the infinite series gives a correction to the 
surface charge density calculated from the previous order 
image. As we discussed after Eq. %M\. the new correc- 
tion induced gets one more factor e~'^^'^' . Since Vs, 
higher order corrections are exponentially small and can 
be completely neglected. 




DR(t>e 



R/r 



(30) 



Let us calculate the interaction energy between a point 
charge and a grounded conducting sphere (Fig. [SJi) for 
the purpose of the next subsection. Considering the in- 
teraction energy between q and q' , we have 



U 



21 e-(n-L')/r'^^-{L~R)/r, 



2D{L - U 
Rq^ 



-2d/r^ 



(31) 



where d = L — R. In the limit of — > c», the exponential 
factor goes to 1 and the energy goes back to the standard 
expression. The additional factor can be under- 

stood in the following way. The charge induced by q on 
the surface of the sphere is proportional to e~'^^^% the 
interaction with it is also proportional to e~'^/^" . 



FIG. 6: The first two image charges induced in each sphere in 
considering the interaction between two spheres. Each image 
charge has its own virtual cloud (not shown) . 



In Sec. mil the energy is calculated using capacitances 
(Eq. (|12|) '). In the present case, it is more convenient 
to consider interaction energy between image charges di- 
rectly. The two methods are equivalent providing that a 
factor 1/2 is added to the interaction energy between 
a charge and its image 7] in the later method (see 
Eq. H36|l l. The two image charges in each sphere are 
given by (see Fig. 



qi = DRcjjie^''' 



DR^h R/r. 



92 = DR^2e'^/'-^ 

, DR^c, 
<l2 = T" 



(32a) 
\ (32b) 



B. Two spheres of the same size 

Now we consider the interaction between two conduct- 
ing spheres of the same size (see Fig. in the presence 
of monovalent salt. Due to the complexity of the method 
of images in the present case, instead of giving the general 
expression for the force, we focus on two limits: c? 3> 
and d <^ Tg. 

Let us first discuss the case when d ^ Vg. As well 
known, when — > oo, there is an infinite series of image 
charges in each sphere joj. For finite r^, similar to what 
we discussed in the last subsection, the magnitudes and 
positions of all image charges are the same but a virtual 
cloud of monovalent salt is created together with each 
of them. In the particular case of d ^ r^, we can actu- 
ally cut off the infinite series and include the contribution 



Here qi and q2 are zero order image charges in each sphere 
which take into account of 0i and (j)2 respectively, similar 
to q" discussed in the last subsection (see Eq. (|30|) and 
Fig. \Sjp) . The screening radius accompanied with them is 
just Ts- The first order image charges q[ and q'2 induced 
by q2 and qi are accompanied by = Rrg/L, similar 
to q' discussed in the last subsection (see Eqs. 
and Fig.lSJt). Similar to Eq. H31I) . the interaction energies 
between these images are 

U{q,.q2) = ^^e-/-, (33) 
C/fe,gi) = -^^e-"^/-, (35) 
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and 
UiL) 



U{quq2) + ^[U{qi,q'2) 



n<P2- 



Correspondingly, 



U{q2.q'i)] 
hi DR^ 



-2d/rs 



(36) 





\ Ts / 




i ' 







F{L) = 010; 



DR^ 



-d/rs 



DR^ 



1 

-2d/rs 



i?2 



L 



L^-R^ 



(37) 



-d/r^ 



appears 
this equa- 



We see clearly that an additional factor of e 
for each order of interaction. When rg — > c 
tion goes back to Eq. H21|l as expected. 

At Ts, the first term in Eq. H37I) dominates. When 
01 =02, we have the standard double layer repulsion of 
DLVO theory P. When 01 and 02 have same signs, only 
the second term represents attraction, which is negligible 
comparing with repulsion. This is analogous to Eq. H21(l 
but attraction is much weaker here. So in the limit of c? 3> 
rs, the force is attractive only if potentials are opposite 
in signs or one of potentials is equal to zero. 

Now let us consider the case when d <C r^. In this 
case, the factor e^^lr^ ^ ^ and all higher order image 
charges should be included. Since each image charge is 
accompanied by its own virtual cloud, the interaction is 
very complicated. Instead of calculating it exactly, we 
estimate C\\ and C12 using a simple method as follows. 
We divide the surface of the two spheres into two pieces. 
In the first piece, the distance between the surfaces of the 
two spheres is smaller than Ts so that the two spheres in- 
teract with each other in an unscreened Coulomb way. 
We call this piece "contact region" (see Fig. Ej). In the 
second piece, the distance between the two surfaces is 
larger than and the interaction between them is expo- 
nentially small and negligible. Therefore, we can safely 
assume that the modification to the capacitances of the 
two spheres due to their proximity happens only in the 
contact region. 

In order to calculate C\\ and C12, let us first consider 
a special case when (f>i — 4>2 — 0- In this case, the 
charge induced on two spheres are equal and given by 
0(C'ii + Ci2)- In the contact region, the screening effect 
by monovalent salt can be ignored and the adjacent spots 
of the two spheres form a capacitor with zero voltage, 
therefore the charge is zero [111. While for the rest part 
of the sphere, the interaction between two spheres is ex- 
ponentially small and the charge on each sphere is almost 
the same as if the other sphere is not there. Therefore 
we have 



Cu + C12 



DR^ _ DRjrs - d) 
rs 4rs 



(38) 



FIG. 7: A schematic illustration of the contact region between 
two spheres. The contact region includes the surfaces of two 
spheres inside a cylinder with the cross section shown by the 
rectangular. 



Here DR^ /rg is the capacitance of a single sphere in the 
presence of monovalent salt. When the sphere gets close 
to the other sphere, it loses its capacitance in the contact 
region, whose area is TTR{rs — d). 

Then we consider another special case when 0i = 
—02 = 0. In this case, we have qi = (Cu — Ci2)0 = — (72- 
Outside the contact region, the capacitance is the same 
as in the previous case. Inside the contact region, the 
capacitance can be estimated as 

dS f^" 2TiR^sinede , rs 

C = D -— = D = 7ri:ii?ln— , 

Jo x{0) Jo d + 2R{l-cosd) d' 

(39) 

where 9o is the angle for the boundary of the contact 
region (see Fig. CJ. The absolute value of the charge in- 
duced on each sphere in the contact region is 2C0. There- 
fore 

Cii-Ci2 = ^-^^(:^ + 2.i.i?ln5. (40) 
Ts 4rc d 



Combining Eqs. (jSHJ and (gUJ), we get 

DR'^ DR{rs - d) 



C12 — 



4r,, 



-7rDi?ln— . 
d 



TTDR\n-j, (41a) 
(41b) 



In the case of rs ~ R, these results match Eq. (|22() which 
was obtained for rs ^ R. 

Using Eq. H14|l , and remembering d ^ rs, we have 



F = -- 



ttDR, 



2d 



h2?. 



(42) 



The force is independent of in this limit because the 
change of capacitances happens mainly in the contact re- 
gion (the logarithm term in Eq. H41|l ). Again, the attrac- 
tion shows up at small d and the attractive force decays 
as a power law. For the model of conductors to work, we 
require d^ A which is valid since ^ <C . 
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C. A small sphere and a big sphere 

Now we consider the interaction between a small 
sphere and a big sphere (see Fig. ^) in the presence 
of monovalent salt. The approach we used is similar to 
that in the last subsection. 

When d ^ Ts, calculating the interaction between 
leading order image charges similar to Eq. (|37|l . we get 

F{L) = M,£^e-'^/^^ 

Again, when and (/>2 have the same signs, attraction 
is exponentially weaker in this limit. When d ^ r^, we 
can estimate Cn, C12 and C22 similar to Eq. H41|l. We 
get an attractive force given by Eq. (Q. 

In the end of this section, let us discuss the possi- 
bility of nonlinear screening by monovalent salt in a 
special case when 01.2 > 0. The result can be generalized 
to other cases. The condition for linear screening reads 

0i,2<0c = ^|ln(c^;o)|, (44) 
e 

where e is the proton charge and c and vq are the con- 
centration and volume of the monovalent ion. Here 0c 
is determined by entropy of monovalent ions in the solu- 
tion. According to Eqs. and (HTJ, 0i,2 ^ Ze/Dri^2 < 
Ze/Da, where a <C ri^2 is the radius of a Z-ion (see 
Fig. I2J). If 01,2 > 0c, Ze/Da ^ 0c. Therefore monova- 
lent counterions are adsorbed to Z-ions. As a result, the 
effective charge of a Z-ion is renormalizcd to 

Z*e ^ Dacfic. (45) 

This effective charge is not changed when a Z-ion is ad- 
sorbed on the surface of the macroion since Z*e/Dri^2 <C 
Z*e/Da. Consequently, 0i,2 should be calculated using 
Z* and satisfy 0i,2 < Z*e/Da = 4>c- In summary, the 
nonlinear screening of monovalent salt is important when 
Z-ions are strongly charged. In this case, our theory is 
still applicable providing that Z is replaced by Z* every- 
where. 



V. ATTRACTION OF TWO MACROIONS, 
WHEN ONLY ONE OF THEM ADSORBS Z-IONS 

In this section, we generalize our theory to a little more 
complicated case when the two bare macroions are oppo- 
sitely charged. As shown in Fig. Positive Z-ions are 
only adsorbed to the negative probe and invert its charge. 
In the presence of large concentration of monovalent salt, 
monovalent counterions are adsorbed to the positive sur- 
face (nonlinear screening), reducing its electric potential 



to 0c (see Eq. (|44ll '). This potential is certainly different 
from the electric potential of the probe determined by 
Z-ions (Eq. dJ). We again have two conductors with 
different fixed potentials. Even when both potentials are 
positive (the net charges of the surface and the probe are 
positive) , the force is still attractive at small distances as 
we discussed in the last section. 

The atomic force experiment in Ref. is actually a 
realization of the situation discussed above. The maxi- 
mum of the repulsive force mentioned in introduction can 
be understood following Eq. H43I) , extrapolated to d ~ Ts . 
When both potentials are positive, at d ^ Tg, the first 
term (repulsion term) in Eq. H43|l dominates and the force 
is repulsive. At d < (r^ ~ 100 A in the experiment), 
the second, attraction term in Eq. H43|) dominates due to 
its larger pre-factor of the exponential. This leads to the 
maximum of the repulsion at d ~ rg. When c? r^, we 
have the polarization attraction given by Eq. (fT|l. 

For d < rs, Eq. ^ can be used to compare the polar- 
ization attraction with the van der Waals one. According 
to discussions in introduction and Sec. we conclude 
that for substantial difference between cti and (72, the 
polarization force Fp is much larger than Fydw at c? 3> A 
(yl = 21 A for Z = 3 and -a = -0.75 e/nm^). Therefore 
the van der Waals attraction has nothing to do with the 
attraction observed in the experiment at d ~ 100 A. 

So far we have considered two cases when the elec- 
tric potentials of the two macroions are different and the 
attraction between likely charged macroions is possible. 
One is that positive Z-ions are adsorbed to both neg- 
ative macroions; the other is that positive Z-ions are 
only adsorbed to one negative macroion, but monovalent 
counterions are adsorbed (nonlinearly screening) to the 
other positive macroion. It is interesting to relate them 
to the more standard case when two likely charged (say, 
negative) macroions are nonlinearly screened by positive 
monovalent counterions (no Z-ions in the solution). In 
this case, the two macroions have the same electric poten- 
tial 0c given by Eq. and never attract each other |^ . 



VI. CONCLUSION 

In this paper, we discussed a long range polarization 
attraction of two likely charged spherical macroions in 
the presence of multivalent counterions (Z-ions). We 
show that the necessary condition for the attraction is 
that the bare charge densities of the two macroions are 
different. This polarization attraction is much stronger 
and longer ranged than the van der Waals force ,5j and 
the short range correlation attraction 6] . In the presence 
of large amount monovalent salt, it adds an additional 
term to the standard double layer repulsion of DLVO the- 
ory when the two macroions are different. We discussed 
two cases when the polarization force between two differ- 
ent macroions can be attractive, even if their net charges 
have the same sign. In the first case, both macroions ad- 
sorb Z-ions (Fig.^). In the second case, one macroion 
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adsorbs Z-ions while the other adsorbs monovalent coun- 
terions (Fig. ^). Here "adsorb" means that the bind- 
ing energy is much larger than fcsT. In both cases, due 
to different equilibrium conditions of adsorbed ions, the 
electric potentials of two macroions arc different and the 
attraction is possible. On the other hand, the attrac- 
tion is impossible if two macroions arc likely charged and 
both adsorb monovalent counterions (no Z-ions in the so- 
lution) . Our result qualitatively agrees with atomic force 
experiments 

Even though in this paper we only discuss spherical 
macroions, the polarization attraction we discovered can 
be generalized to other geometries. Actually, as seen be- 
fore, the attraction is always developed at small distances 
between macroions when the overall geometry is not very 
important. One can also understand this using the lan- 
guage of contact electrification as discussed in the intro- 
duction. 

In this paper, in the connection with atomic force ex- 
periments [2, 13 J we focused on the force between two 
macroions instead of the total free energy of the system. 
Actually, for typical situations (e.g., A <^ R oi A <^ Vs), 
one can show that the free energy may have the global 
minimum when the two macroions are close to each other. 
Thus, the polarization attraction can play an important 
role in determining the equilibrium state of the system. 
It answers the gene delivery related question: even if 
charges of DNA and cell membrane are both inverted, 
they can still be attracted to each other since DNA and 
membrane certainly have different surface charge densi- 
ties. It can also be important for aggregation or self- 
assembly of large ensembles of different likely charged 
macroions with help of oppositely charged Z-ions. As- 
sume for example that we have a mixture of equal num- 
bers of two kinds of negative spheres with the same radius 
but different values of the bare charges. Then in the pres- 
ence of positive Z-ions, they can attract each other and 
assemble in a NaCl-like structure. 



Evaluating the integral, we have 

1 , cosh t + 1 ■r^ 1 1 , cosh t + 1 1 

— In <> <— In \ . 

2t cosh t — 1 ^-^ sinh nt 2t cosh t — 1 sinh t 

n—l 

In the limit of i ^ 0, it becomes 



(A.2) 



- In - < > — < - In — . (A.3) 

t t ^ sinhnt t t ^ ' 



1 1 . 2e 



Therefore in this limit, we have 



sinh nt t t 



(A.4) 



where /? is a number satisfying 1 < /3 < e. Replacing t 
by 2t, we get 



°° 1 



iln^ 
sinh 2rd 2t t 



(A.5) 



And 



y - = y -J: y t 

sinh(2n — l)t sinhnf sinh 



n=l 



2nt 2t t 
(A.6) 

Since coshi — {2B. + d)/2R, the limit of t ^ is equiv- 
alent to d ^ and t — y/d/R in this limit. According to 
Eq. l(T3|l . we have 



DR^ WP'^R 
DR P^R 



(A.7a) 
(A.7b) 



APPENDIX 



Numerical calculation gives /3 — 1.78. 



Let us start from the asymptotic behavior of the se- 
ries ^2"^=! 1/ sinhnt at t ^ 0. According to the Cauchy 
Integral test of convergence, we have 
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